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We analytically and numerically analyze the one-dimensional "thin-torus" limit of Fractional 
Topological Insulators in a series of simple models exhibiting exactly fiat bands with local hopping. 
These models are the one-dimensional limit of two dimensional Chern Insulators, and the Hubbard- 
type interactions projected into their lowest band take particularly simple forms. By exactly solving 
the many-body interacting spectrum of these models, we show that, just like in the Fractional 
Quantum Hall effect, the zero modes of the thin-torus limit are CDW states of occupation numbers 
satisfying generalized Pauli principles. As opposed to the FQH where the thin-torus CDW appear in 
orbital space, in the thin-torus FCI states, the CDW states are in real-space. We show the counting 
of the quasihole excitations in the energy spectrum cannot distinguish between a CDW state and 
a FQH state. However, by exactly computing the entanglement spectrum for the thin-torus states, 
we show that it can qualitatively and quantitatively distinguish between a CDW and a fractional 
topological state such as the FCI. We then discover a previously unknown separation of energy 
scales of the full FQH energy spectrum in the thin torus limit and find that Chern insulator models 
exhibiting strong isotropic FCI states have a similar structure in their thin-torus limit spectrum. We 
close by numerically computing the evolution of energy and entanglement spectra from the thin-torus 
to the isotropic limit. Our results can also be interpreted as an analysis of one-body, 1-dimensional 
topological insulators stabilized by inversion symmetry in the presence of interactions. 
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I. INTRODUCTION 

Nontrivial topological insulators are systems not adia- 
batically continuable to an atomic limit through a path 
obeying some set of symmetries (such as time-reversal, 
charge conjugation, or a point group symmetry) which 
characterize the topological clas^l. When an edge or 
surface is present, the topological insulator may or may 
not exhibit gapless edge or surface modes depending on 
whether the symmetry characterizing its bulk is unbroken 
or broken by the presence of the surface. The theoreti- 
cally proposecPlSI and experimentally discoverecP^ topo- 
logical insulators are time-reversal invariant, and can be 
explained mainly through one-body theory, as their par- 
ent materials exhibit weak interactions. However, recent 
theoretical research has focused on proving the existence 
of new interacting states of matter, the Fractional Topo- 
logical Insulators (FTI), which come about when one or 
more of the topologically nontrivial bands of a one-body 
nontrivial topological insulator have narrow bandwidth, 
are fractionally filled, and are subject to interactionJ^tHl 
When the interactions overcome the weak band disper- 
sion, new topological phases occur. The parade example 
of an FTI is the Fractional Chern Insulator (FCI) : it 
consists of interacting electrons in a fractionally filled, 
specially designed fiat band (zero bandwidttP^^^^^ with 
nonzero Chern number^''. For some specific values of the 
filling factor, they form a collective state s imilar to that 
in the Fractional Quantum Hall effeclP^^ but in the ab- 
sence of an overall magnetic field (fiux) on the sample. 

Although several attempts have been made to under- 
stand the occurrence of FTI states from first principles, 
many unresolved issues remain, even in the simplest ex- 



ample of an FTI, that of the FCI. First, it is unclear 
how to differentiate the FCI state from broken symme- 
try states such as charge density waves (CDW) with small 
order parameter: on the lattice, a commensurate CDW 
state is gapped, could appear at identical momenta as the 
FQH state, and could be favored by Umklapp processes. 
The arguments based on the counting of the quasihole ex- 
citations in the energy spectrum in the FCI put forward 
by could also be applied to the CDW state. It is also 
unclear which interactions favor a FQH state in specific 
one-body models: while a correlation between s mooth 
Berry curvature and the existence of the FCI stat^^lllll 
has been observecP^, we still do not know how (and in 
which models) the interplay of local interactions in real 
space with the local eigenstate structure in momentum 
space leads to robust FCI states. 

In the FQH effect on the toroidal geometry, useful in- 
formation can be extracted from a limiting procedure 
called the "thin torus limit"^^^ in which the aspect ra- 
tio of the torus is skewed such that one length \Li \ of the 
torus is taken to zero while the product Li ■ L2 — 2nN^ 
is kept constant, with equal the number of fluxes 
through the system. In t his lim it, convincing analyticaP^l 
and numerical argumentJpfil2l! have shown that the FQH 
ground-states and quasihole excitations are just the sin- 
gle slater determinant of orbital occupation numbers ex- 
hibiting the possible patterns of commensurate CDW at 
the filling N^/N^ = p/q with p, q relatively prime and 
short-range Hamiltonians given by electro-static inter- 
action (in the case of the Laughlin states) or a gen- 
eralized fc-body electrostatic interaction in the case of 
the Zk Read-Rezayi (RR) statePH^. These CDW pat- 
terns are on the torus are identical to the "root parti- 
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tions" pattern^^H which index the many-body isotropic 
FQH states. As the torus is deformed smoothly from 
the \Li\ to the \Li\ = IZ2I isotropic Hmit, it was 
showed that the CDW states evolve continuously into the 
isotropic FQH states which exhibit no CDW order pa- 
rameter. However, the excitations of the isotropic FQH 
states can be counted in the same way as the excitations 
of the CDW thin-torus states, even though they exhibit 
no local order parameter. . 

In this paper we perform a similar "thin-torus" anal- 
ysis of the Fractional Chern Insulators by analyzing the 
FCI problem on a lattice with — 1 sites in the x di- 
rection. While on the lattice it is not possible to continu- 
ously interpolate between the thin torus and the isotropic 
limit, subsequent increases of the from 1 to 2,3,... 
while keeping N^Ny constant (when possible) show how 
the thin torus CDW states evolve into the FCI states in 
the isotropic limit. We propose a series of simple lat- 
tice 1-D models which have an exact flat band with local 
short-range hoppings. These models are the 1-D limit 
of 2-D one-body Chern topological insulators, to which 
we then add Hubbard-type interactions. The full spec- 
trum of these interacting one-dimensional models in the 
topologically trivial and nontrivial sides can be obtained 
exactly for any filling in both the topological side and 
the atomic limit. The quasihole spectrum, which in this 
case is made up of zero modes, exhibits the counting of 
excitations of the FQH state, even though the state is 
a CDW, and hence cannot be conclusively used to de- 
termine whether an FCI state in the isotropic limit is in 
the same universality class as the FQH. We however find 
that the structure of the full energy spectrum in the thin 
torus limit, including excitations, exhibits a separation of 
scales in the FQH effect into different bands of states dis- 
tinguishable by their electrostatic energies. We observe 
the same separation of scales in the thin-torus limit of 
Chern insulator models with strong isotropic FCI states. 
The particle entanglement spectrunPH of the (degener- 
ate) ground-states can also be obtained exactly in the 
1-D limit. We show that the entanglement spectrum dif- 
fers considerably, both quantitatively and qualitatively, 
in the thin torus CDW limit from the isotropic FCI limit, 
and hence can be used as a diagnostic for the topological 
nature of the state. We then numerically analyze several 
different models which do not have a flat band in the thin 
torus limit, and show that the structure of their spec- 
trum evolution from the thin-torus to the isotropic limit 
matches that of the FQH evolution between the same 
limits. Our results show that a 1-D one-body topolog- 
ical insulator stabilized by inversion symmetry exhibits 
CDW states upon the introduction of local interactions, 
and that the only difference between the trivial and non- 
trivial Inversion symmetric interacting topological insu- 
lators is the range of the interactions when projected to 
the lower band. 

The paper is organized as follows. In section 2 we de- 
scribe the general framework for the one-body models, 
including the flat band limit. We also give a detailed of 



the two orbital and Kagome models. Section 3 is devoted 
to the interacting two orbital model. We show that its 
one dimensional limit can be handled analytically in the 
isotropic limit. We prove that at filling factor ly — 1/3, 
the low energy physics is a CDW. In Section 4, we com- 
pute that the particle entanglement spectra of a CDW 
and we show that it differs dramatically from the one of 
the Laughlin state. Section 5 deals with the energy spec- 
trum structure of the FQH in the thin torus limit. We 
numerically study several FCI models in the thin torus 
limit and compare them to the FQH case. 

II. THE ONE-BODY MODELS 

In this paper, a one-body translational invariant A^- 
band Bloch Hamiltonian H ~ j cl^hap{i—i)cji3 where 
the i and j denotes the unit cell coordinates, a, P = 
1,...,N are orbital (and spin, when appropriate) in- 
dices, the summation over a and (3 is implicit. The 
Bloch Hamiltonian is named local if its matrix elements 
ha/sO- ^ j) decay exponentially with the distance i — j. 
Usual tight-binding Hamiltonians are of this form, with 
hopping taking place usually between close-neighboring 
sites. The Fourrier transform H = J2]i'^a^ap0^)c\ij3 
of such a Hamiltonian contains only limited harmonics 
of the Bloch momentum k. The model can be diago- 
nalized in normal modes 7"(k) of energy En(k.) where 
n — 1, . . . ,N is the number of energy levels at each 
i'-dimensional (with D — 1,2,3 experimentally acces- 
sible) momentum k = {2TTii/Ni, . . . ,27Tio/No) where 
ij = 0, . . . Nj — 1 runs through the number of sites in 
the j — 1, . . . , Z3'th direction. The diagonalized Hamilto- 
nian reads H = X^n k -^n(^)7"^(^)T"('*) where 7"(k) — 
u^^Ckjci^a and u^(k) diagonalizes the Bloch Hamiltonian 
matrix hap(k)u^{k) — i?„(k)u2(k). The projector oper- 
ator in the n'th band reads = 7"1'(k) |0) (0|7"(k) 
with matrix elements Pj^^ (k) = u2(k)u^^(k). In terms 
of projectors, the one-body Hamiltonian can be written 

as/i„Mk)=E„^n(k)i;'i?(k). 

If the one-body Hamiltonian analyzed is an insulator 
in which a number Ni of the "low" energy bands are 
separated by a full (direct and indirect) gap from the 
rest " high" energy N — Ni bands, we can then introduce 
a deformation H^^ of the Hamiltonian H which keeps 
its eigenstates intact but changes its energies so that they 
are non-dispersive or Flat Band: 

- E E -^^'"(k) + E E^^"Hk) (1) 

n=l k n=Ni k 

The simplified fiat-band Hamiltonian and not the full 
Hamiltonian is sufficient to analyze the topological 
properties of an insulator. For example the Chern 
number of the occupied bands can be expressed as 
C = EkeyTr[PG(k)a,PG(k)a,PG(k) where PG(k) = 
X^j^j^ p("^(k) is the projector into all of the occupied 
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bands. The simplified flat-band Hamiltonian is now 
equivalent in spirit to a Landau Level in the sense that its 
zero kinetic energy does not compete with with any inter- 
action. At fractional filling of any flat band, the one-body 
system is a zero-bandwidth metal and exhibits thermo- 
dynamic limit degeneracy of the metallic states. Upon in- 
troducing repulsive interactions which are projected into 
the fractionally filled band (thereby artificially render- 
ing gaps infinite - a limit similar to that of projection in 
the Lowest Landau Level) , two things can happen. In 
the thermodynamic limit the degeneracy is lifted, but, 
supposedly in different ways depending on whether the 
one-body system is topologically trivial or nontrivial. In 
the atomic limit, a repulsive interaction will still (most 
likely) allow for a large degeneracy of states characterized 
by the configurations of particles on sites minimizing the 
repulsive local energy. On the side where the one-body 
Hamiltonian is topological nontrivial, the hope is that the 
delocalized nature of the eigenstates will frustrate the lo- 
cal interaction such that the system will exhibit a finite 
degeneracy in the thermodynamic limit. A further hope 
is that this degeneracy is then not due to a symmetry 
broken gapped state on the lattice such as a CDW, in 
the thermodynamic and isotropic limit. 

Several Chern insulator models have been studied in 
presence of st rong repulsions : the checkerboard lattice 
modelJ ^^° l ^^ l, the Halda ne honeycomb model^ iHM I the 
ruby lattice mode^^^, the Kagome lattice modeP^'^^ 
and the two orbital model on a square Mtic^. In this 
article, we will mostly focus on these two latest examples 
that we will now describe. 

The simplest example of a Chern insulator model is 
a two-orbital model on a square lattice with s and p 
orbitals on site as depicted in Fig. [T^. The nearest 
neighbor Hamiltonian with a topological transition reads 
H = Ek(4,s'Ck,p)'*s-p(k)(ck.s,Ck,p)^ : 



/is_p(k) ^ ^ -sm{ki)(Ji - (Af 



t—x,y 



i—x,y 



cos{ki))a2 (2) 



where (Jx,y,z are the three Pauli matrices. This Hamilto- 
nian has 4 phases, two adiabatically continuable to the 
atomic limits M — ±oo with Chern number C — Q, one 
phase < Af < 2 with Chern number C = — 1 and one 
phase — 2 < M < with Chern number C = 1. M = 
2, 0, —2 are one-body topological phase transition points 
where the gap closes in a Dirac fermion mass-changing 
transition at k = (0, 0), ((0, tt) and (tt, 0)), (tt, tt). In fact, 
the above Hamiltonian Eq[2] is the first of a part of a 
series of Hamiltonians /i"_p with n'th neighbor hopping: 

/i"_p(k) = sin{nh)a^ + {M - ^ cos(nfc,))cr^ (3) 

i—x,y i—x,y 

which also exhibit 4 phases, two adiabatically continu- 
able to the atomic limits M — > ±oo with Chern num- 
ber C — phase < M < 2 with Chern num- 



ber C 



and one phase — 2 < Af < with 



topological phase transition points where the gap closes 
in a Dirac fermion mass-changing transition at k = 
{2TTii/n, 2ni2/n2), ((27rii/n, 2n{i2 + l/2)/n) and (27r(ii-|- 
l/2)/n, 2m2/n)), {2n{ii + l/2)/n, 2n{i2 + l/2)/n) where 
«i, ■*2 = 0, . . . , n — 1. 





FIG. 1: Left panel: The two orbital model on the square 
lattice. The two orbital s and p in each unit cell are depicted 
by the gold and silver sphere. Left panel: the Kagome lattice 
model made of its the three sublattices A, B, C. The arrows 
give the sign convention for the phase of the NN hopping 
term; the hopping amplitude is exp(iiy9) in the direction of 
the arrows. The lattice translation vectors are ai and 02. 

The Hamiltonian Q have a discrete point group sym- 
metry under inversion V = cJz'- 

Vhl^{\.)V-'=hl^{-\,) (4) 

The inversion eigenvalues C(K) of the lower 
band at the inversion symmetric points K = 
(0, 0), (0, tt), (tt, 0), (tt, tt) change in the phase transition 
change in a phase transition, and the Chern number par- 
ity is consistent with the relation (—1)*^ — Hk C(K)r^. 

The other model that we will focus on, especially for 
the numerical calculations, is the Kagome modeP^l. It 
consists of three sublattices A, B, C spanned by the 
translation vectors ai and 02 (see Fig [T]d). Different 
tight binding models leading to a nonzero Chern number 
have been consider in Ref.E^. We consider the simplest 
case with only the nearest neighbors hopping terms of 
amplitude exp{±iip). The Bloch Hamiltonian reads 



h{k) 




h.c. 



(1 



-iki 



'(1 



(1 



Chern number C 



M 



2, 0, —2 are one-body 



(5) 

where fci ~ k ai and ^2 = fc-£i2- The phase diagram of 
this model is determined by the (f phase. The gap closes 
at If = 0, ±^ and away from these points the lowest 
and highest energy bands have a Chern number equal to 
±1. When looking at the interacting case at filling factor 
V — 1/3, this system has been shown to host the most 
robust fermionic Laughlin-like stat^^^. So it is the most 
natural candidate to consider if one wants to address the 
question of the thin torus limit. As we will see below, its 
spectrum in the thin-torus limit is very similar to that of 
the FQH state without the need for fine-tuning. 
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III. THE INTERACTING TWO ORBITAL 
MODEL 

We now add a density-density interaction ranging from 
on-site to the m'th nearest neighbor: 

= ■ ('^'^^ + '^i,p)("j,s + : (6) 

ij;|i-j|<™ 

where ni^s,ni_p are the densities Cj ^Ci^s, Cj ^Q^p of the on- 
site s,p orbitals and where the term : . . . : denotes the 
normal ordering usually performed in FQH calculationJ^Sl 
where all the f operators get placed to the left with only 
signs due to permutations. We restrict to V|i_j| all pos- 
itive. Note that for the orbital model this interaction 
contains an on-site nj.srij.p interaction between the two 
orbitals. The interacting Hamiltonian is called local if m 
is finite with m/Nx — >■ as — > cx). To work in the low- 
est band, we now project this interaction to the negative 
energy band (this corresponds to setting the interaction 
energy much lower that the artificially modified one-body 
gap, which is now infinity). In momentum space the pro- 
jection is local at every k (as the one-body Hamiltonian 
is diagonalized in k). In real space, however, the interac- 
tion now becomes long-ranged (the interaction will have 
matrix elements between any two sites on the lattice no 
matter what the number of sites , Ny is) . If the system 
is gapped, the interaction matrix elements will still de- 
cay exponentially for |i — j| > m, where m is the original 
range of the Hamiltonian. Projection of the Hamiltonian 
to the lowest band has rendered the problem difficult as 
the eigenstates are strongly dependent on k, especially 
in the topological nontrivial side, where the "winding" of 
the eigenstates gives rise to the Chern number. 

A. Isotropic Limit: Overview of General Properties 

One limit of the 2D problem is easy to analyze: in the 
atomic limit A/ — ^ cx) (the limit M — oo is of course 
identical) the projection to the lowest band fixes the 
problem to contain only s orbitals (the p-orbitals have an 
infinite energy). The projected interacting Hamiltonian 
becomes P^^^Hi^P^^^ = Ei,j;|i-j|<m : ni.sjij^s ■= 

^i7^j |i-j|<m ^i^l ■ '^i-s^'j.s •• spcctrum Can be deter- 
mined by combinatoric counting. For example for on-site 
interaction m = in the unprojected Hamiltonian, the 
projected Hamiltonian becomes zero and the spectrum 

consists of ^ ^xNy ^ ^^^^ ruodes which count the con- 
figurations of Ne fermionic particles in N^Ny boxes. 

We now ask what happens in the interacting prob- 
lem when the one-body Hamiltonian is driven through 
a phase transition from the atomic limit M -> oo to the 
< M < 2 phase with Chern number nonzero. With 
interaction projected into the lowest band, numerically, 
at least for A'^j, x A'^y = 6 x 4 and 6 x 5 at iVg = 8, 10 
electrons respectively (fiUing 1/3), a 3-fold degenerate 



ground-state is observecP^ with on-site Vq and nearest 
neighbor Vi interactions, but only for Vi « Vq. For a 
yet undetermined reason, no gapped degenerate ground- 
states is obtained for smaller sizes. The observed states 
are fragile with respect to Vi (see Figs. [2^ and b), though 
a ratio Vi/Vq = 0.0015 leads to a 40% increase of the gap 
over the Vi = value. The fact that the optimal Vi/Vo 
ratio is small will become important later. 

B. Thin Torus Limit: One Body Model 

While we cannot solve the interacting model analyti- 
cally in the generic situation, the thin torus limit — 1 
is amenable to theoretical analysis. In this limit the 
Hamiltonian becomes: 

/i"_p i^(k) = — sin(nfc)a-j, — {M — 1 — cos{nk))(Tz (7) 

For n odd and < M < 2, this Hamiltonian is a topo- 
logical insulator stabilized by inversion symmetry with 
P = az- For AI > 2, the low energy band has inversion 
eigenvalue ^(0) = CC'"') = ~1 for all n (odd and even). 
At A/ = 2, a phase transition flips the inversion eigen- 
values at A; = only for n odd but at both k — 0,tt 
for n even, and these eigenvalues remain identical in the 
gapped phase until M — when the inversion eigenvalue 
at A; = TT flips for n odd, while for n even there is no in- 
version eigenvalue change as the gap closing (for AI — 0) 
does not take place at an inversion symmetric point in 
the Brillouin zone. Hence in the case of n odd, the prod- 
uct of eigenvalues C(0)C(^) = — 1 for < M < 2 and the 
system exhibits a charge polarization J dkAi^/ii: = 1/2 
( mod 1) where = —iu^a'^^'^ka Berry phase 

of the lower band with eigenstate . The value of the 
charge polarization is fixed by inversion symmetry and is 
an invariant of the one-body topological phase - it cannot 
be changed unless a phase transition occurs. By contrast, 
the n even models exhibit unit charge polarization and 
are trivial under the inversion symmetry charge polariza- 
tion classification. When the system is cut in half, the cut 
(or edge) breaks inversion symmetry, and the system does 
not need to exhibit mid-gap modes in its open-boundary 
energy spectrum (although the Hamiltonian ^ does ex- 
hibit these modes, one on each cut, due to the absence of 
a kinetic term e(fc) proportional to the identity matrix, 
which we neglected for aesthetic purposes but which is 
allowed by inversion symmetry). The system however 
exhibits mid-gap states in its entanglement spectrunP3^ 
another way in which it reveals its topological nature. 
For the 2-dimensional models we chose, their thin torus 
limit is, for n odd a topological insulator with inversion 
symmetry, and our analysis also reveals the fate of an 
inversion symmetric topological insulator in 1 dimension 
(we note that one needs a symmetry in ID to define a 
topological insulator). Our results, however, will be valid 
when inversion symmetry is absent, and we focus on the 
thin-torus limit of the FCI rather than on the 1-D inter- 
acting inversion symmetric topological insulator. 
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We note that the Hamihonian ([8]) makes sense only if 
Ny/n S Z. By contrast, the trivial atomic limit Hamil- 
tonian reads: 



H 



M- 



(10) 



Note that the expression of the Hamiltonian in terms 
of the new fermionic operators aj,&j is just the rotated 
Hamiltonian Eq[7]: to the sin(nfc)crj, + cos{nk)(Jx basis. 
A simple decomposition into normal modes is also more 
revealing in real space: 



V2 



V2 



(11) 



,t _ uf 



0-j bj+n flj - hj+n 



72 



72 



and 



FIG. 2; Upper panel : Energy spectrum for N — 8 fermions on 
the two orbital model with Ny — 6, Nx = 4 and Af = 1. Only 
the two lowest energies per momentum sector are displayed. 
Both spectra for Vi and Vi/Vo = 0.01 are shifted with their 
respective ground state energy Ei. For this system size and 
this aspect ratio, the Laughlin state v = 1/3 should appear in 
momentum sectors {Kx,Ky) = (0,0), (0,2) and (0,4). Lower 
panel : Gap between the fourth state energy and the third 
state energy as a function of the nearest neighbor repulsion 
Vi. Until the gap closes, the three lowest energy states have 
the momenta corresponding to the Laughlin state. 



We further note a special point of the model in Eq[7]: 
the point for which the gap is maximal is also the point 
at which the Hamiltonian has an exact flat-band without 
the need for projection. Indeed, for Af = 1, the energy 
spectrum of the Hamiltonian is E± = ±1. Note that 
only the thin torus limit has a fully flat band at Af = 1 - 
the isotropic limit has dispersion in the second direction. 
Generically, however, it is true that the thin torus limit 
Hamiltonian has flatter bands than its isotropic coun- 
terpart. The fact that the Hamiltonian /i"_p]^^(k) in 
Eq[7] has a flat band with only short range (n'th neigh- 
bor maximum) hopping is revealing as it means that the 
projectors into the lowest energy band are short ranged 
and hence amenable to analytic derivation. The M = 1 
Hamiltonian is easiest to understand in real space: 



H 



- Ek(Ck,s' Ck,p)(sin(nA:)(Tj, + cos(n/c)cr^)(ck,s, Ck,p)'^ = 
= - Ejii a]bj+„ + h.c. (8) 

where the new fermionic operators aj , bj on site j are a 
combination of the s,p on-site orbitals: 



1 



[Cjs Cjp), 



1 

71 



{Cjs + Cjp) 



(9) 



M- 



- 2^7 = 1 



'/2 \/2 



LAiLiL (12) 



\/2 \/2 



The normal mode of the lowest band is hence 



7j 



{ttj + 6j_|_„)/%/2 on the nontrivial and 7^ = {aj + bj)/\/2 
on the trivial side respectively. Hence the nontrivial side 
is a generalized Su- Schrieffer-Heeger model with perfect 
n + 1-merization. The ground-states of the two limits 

are V'm=i ~ Il^il'^j ~^ ^j+n)\0) containing dimeriza- 
tions between site j and site j + n on the nontrivial site 

and i^M^oo — Ylf=ii^j + ^j) |0) s-wave orbital 

on each site j occupied. The most well-known model of 
the non-trivial side is the n = 1 model which contains 
a dimerization between site j and j + I, with the center 
of charge mid-bond. A cut on the bond cuts the charge 
center. The n-generalization of this model exhibits ex- 
pectation values (dimerization) of operators on site j and 
j + n. Cutting one bond cuts through n of the dimerized 
lines. 



C. Thin Torus Limit: Projected Interaction 

We now try to project the interacting Hamiltonian 
to the lowest bands. In general, in the isotropic limit, 
the projection to the lowest band has several conse- 
quences. First, even though the interacting Hamiltonian 
Eq[6] is density-density, projection to the lowest band 
will make it unwritable in terms of just density oper- 
ators on site. Second, it will introduce nonlocal long- 
range couplings that run through the length of the lat- 
tice. This makes the Hamiltonian very difficult to ana- 
lytically investigate. The current 1-D flat band case is 
different. Due to the fact that the projection is local, the 
projected Hamiltonian is also local. The projection to 
the lowest band of any interacting Hamiltonian is now 

easy to do and accomplished by taking Pn^'ajPn'' — 

7j/-\/2, Pn^bjP!^^ = ^j-nl\fi in the nontrivial phase 
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and Pi^^ajP^P = 7^/72, p[^h,Pi^^ = jj/V2 on the 

trivial side, where Pn^'' is the projector into the lower 
band for the Hamiltonian /i"p]^^(k) at the points AI — 1, 
M — >■ oo (nontrivial side and atomic limit). Hence the 
projected initial operators of the s and p orbitals in the 
nontrivial phase read P^^ Cj^P^^ — (7j + 7j_„)/-\/2 and 
P^i'^ CjpP^^ — {'Jj-n — 7j)/\/2, while in the atomic limit 

they are Pi^'^c^sPn^ = Ij and Pk^^ CjpPi^^ = 0. A 
density-density term interacting between r sites in the 
original Hamiltonian Eq[6] is nj^ra,nji3. When normal 
ordered and projected to the lowest band becomes, we 
obtain for n ^ r : 

= ^ 2nj+rnj + fljj^r+nflj + hj+r-nflj (13) 
3 

for n = r > : 

: ^ P^'rij+r.anj^pP'^ : 
= ^ 2hj+rnj + rij+2r%- (14) 

3 

and for r — 0: 

where fij = 'Jjjj is now the density operator of the nor- 
mal modes jj. One sees that in the topologically or- 
dered side, an interaction between the particle densities 
r sites away becomes, when projected to the flat band of 
/i"pj£,(k), an interaction at r -I- n sites away. The only 
result of the projection is to increase the range of the in- 
teraction. By contrast, on the trivial side, the projected 
Hamiltonian is zero for r = and for r 7^ - 

in other words, the range of the Hamiltonian remains 
unchanged. The density-density projected Hamiltoni- 
ans act on the Hilbert space Hil^^i 7ji and solving 
them becomes as simple problem of solving for config- 
urations minimizing a hard-core density-density interac- 
tion with no kinetic energy. The interactions can now be 
tuned to obtain the correct degeneracy of FCI states in 
the thin torus limit, which is known from the FQH thin 
torus limit. In that case, the Hamiltonians are density 
-density in momentum space (as we will see in the next 
sections) and its degenerate ground-states are zero-mode 
single Slater determinants of CDW-patterns in momen- 
tum space^'^^^^ . It is now clear that all the FCI thin torus 
states are commensurate CDW states but in real space. 
A similar analysis can obviously be extended to higher- 
density interactions such as those necessary to give a 
Moore- ReacP^ or RBpH state. For bands with some dis- 
persion, which are not flat-band from the get-go, other 
longer range terms of the Hamiltonian are generated, but 
we expect that the major properties of the spectrum are 
maintained. 



D. The = 1/3 case 

We now particularize to the case of electron filling 1 /3 
and Chern number 1 in the isotropic limit. This means 
we are looking at the /i"p Hamiltonian which for M = 1 
and Nx = 1 has the thin torus limit hl_^^jj with a flat 
band. To obtain a 3-fold exact degeneracy (same as that 
of the Laughlin FQH state) in the thin torus limit it 
is now obvious one needs a combination of on-site (inter- 
orbital) and nearest neighbor density-density interaction: 

Hi/3 = ^0 ^ ?^i,s'^i,p + Vi ^(^i^^ + ni^p){nj^s + n-j,p) 

i <iJ) 

(16) 

where normal ordering is no longer an issue. The relative 
scales Vi/Vq for which a FQH state is seen in the isotropic 
limit will be borrowed from the analysis of the FQH prob- 
lem in the thin torus limit, to be presented soon. When 
projected to the lowest band, in the topologically non- 
trivial flat-band case M = 1, this Hamiltonian becomes: 

P^^^ Hi/aP^^^ = {Vo + 2Vi)J2f^ini + Vi ^ ninj (17) 

(iJ> ((iJ)) 

while in the atomic limit M ^ oo the Hamiltonian be- 
comes: 

=^i^mnj (18) 

(iJ> 

Notice that the difference between the projected Hamil- 
tonians in the non-trivial and trivial phases is the range 
of the interaction. At a certain filling, the non-trivial 
Hamiltonian will have far less zero modes than the triv- 
ial one. 

The spectrum of[T7|is easily obtained for Vo,Vi > 0. 
At filling Ne/Ny = 1/3 {Ny is the number of sites in 
the ID lattice), the Hamiltonian has a 3-fold degenerate 
groundstate which is a charge density wave 

Afe A^e 

|1) - n^3. |0) , |2) = n^s^+i |0) , |3) =117^+2 |0) , 

i—1 i—1 i—1 

(19) 

These are the thin-torus CDW states in which the Laugh- 
lin state evolves. At lower filling N^/Ny < 1/3 there are 
many zero modes given by the number of ways of putting 
Ne electrons in Ny sites without any two electrons com- 
ing closer than 3-sites together (in the FQH literature, as 
mentioned below, this is the (1,3) Pauli principle, though 
there it appears in orbital space). 

The excitation spectrum of [Tt] can also be easily ob- 
tained exactly, although for the case when Vq, Vi are of 
the same order of magnitude, it involves a certain amount 
of book-keeping. However if, as in the FQH thin torus 
case analyzed in the next section, Vq >> Vi (what we 
need is that Vq >> {N^ — l)Vi, which is more than sat- 
isfied in the thin-torus FQH), the excitation spectrum 
separates into two bands, a phenomenon which we call 
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separation of scales: one band, between energy Vi and 
{Ne — l)Vi is the band of wavefunctions which are Slater 
determinants 0^=1 1^) such that the distance between 
any two sites 3 > | j^^ — ji^ \ > 1 and at least one pair 
of j's have distance 2 (i.e. not all j's be separated by 
a distance 3, as those would be zero modes). In the 
FQH, these configurations of Slater determinants satisfy 
a (1,2) Pauli principle which says that there can be no 
more than 1 particle in 2 consecutive sites. Notice that 
the (1,3) Pauli principle is a sub-principle of the (1,2) 
Pauli principle, and hence the number of excitations with 
energies between Vi and {N,, — l)Vi is equal to the num- 
ber of configurations of electrons in Ny sites satisfying 
a (1,2) Pauli principle minus the number configurations 
satisfying a (1,3) Pauli principle. The last band of ex- 
citations, from energy Vq + 2Vi to {N^ — l)(Vo -I- 2Vi) is 
that of all fermionic configurations (Pauli principle (1, 1)) 
minus the (1,2) and (1,3) configurations. This is the 
nontrivial M = 1 case. The current explanation of the 
spectrum can be extended to other Hamiltonians such as 
those that give the RR states - with similar conclusions 
- in the thin torus limit the zero modes are given by the 
CDW patternPSMU or the root partition patterns^^' of 
the respective RR state in real space. The atomic limit 
M — >■ oo has a spectrum made out of all 1,2 excitations 
as zero modes (it has a size-dependent and increasing 
number of zero-modes) and the remaining (1,1) excita- 
tions at energies Vi to {Ne — l)Vi (not Vq). Note that our 
analysis of the current 1-D system can also be interpreted 
as the analysis of a topological insulator stabilized by in- 
version symmetry, although the ground-states we obtain 
are always charge density waves. 



The thin torus limit of this exactly solvable model 
teaches us several things: first, the PCI states in this 
case are exactly degenerate zero modes trivial commen- 
surate CDW state in real space (3 states for Laughlin 1 /3 
state). Second, the number of zero mode excitations at 
lower filling than 1 /3 (the number of " quasiholes" ) is the 
same as in the isotropic FQH, even though in the thin 
torus limit they are excitations of the trivial CDW state 
(they are single Slater determinants). This last point 
shows that it is not possible to infer the character of 
a FQH state just by counting excitations in the energy 
spectrum - a CDW state would also have the same ex- 
citations. Third, in order to mimic the behavior of the 
FQH spectrum (obtained below), the interaction has to 
be such that Vq >> {N^ — l)Vi. Indeed, in the two- 
orbital model [2] turns out to only have a 3-fold degen- 
erate state separated by a gap in the isotropic limit for 
Vq >>> {Ne — l)Vi as we have seen in the numerical 
data in Figs. [2^ and[2]D. In this case, the thin torus limit 
spectrum of the FCI and FQH is very similar (including 
the separation of scales mentioned above and in the FQH 
section below). 



IV. ENTANGLEMENT SPECTRUM 

The counting of the quasihole subspace cannot be used 
as a clear diagnostic of the FCI state in the isotropic limit, 
as the CDW state (which is actually the ground state 
in the thin-torus limit) has excitations obeying identi- 
cal counting to that of the FQH quasiholes. How can 
we then clearly identify that the state in the isotropic 
limit is a topological (FQH) state and not a CDW (with- 
out computing order parameters)? We now show, by 
direct computation, that the entanglement spectrurrp2l 
can distinguish between the CDW state and the isotropic 
FQH/FCI topological state. We use the particle en- 
tanglement spectrum^'^ (PES). For a d-fold degenerate 
state {{ipi >}, we consider the density matrix p — 

1 J2i=i IV'i X V'il- We divide the N particles into two 
groups A and B with respectively Na and Nb parti- 
cles. Tracing out on the particles that belong to B, we 
compute the reduced density matrix pa = Tig p. This 
operation preserves the geometrical symmetries of the 
original state, so we can label the eigenvalues exp(— ^) of 
PA by their corresponding momenta, but generically we 
can also look at all momenta sectors without separating 
quantum numbers (which is what we will do here). The 
entanglement spectrum is just the ^'s (generally called 
energies) plotted as a function of the momentum. Gener- 
ically, for a contact (such as pseudopotential) interaction, 
the number of non zero eigenvalues in pA is bounded from 
above by the number of zero-modes of the Na particles 
in the large Hilbert space of the original state. For ex- 
ample, for a pseudopotential interaction in the FQH ef- 
fect for Ne electrons in N^ fluxes, the maximum number 
of nonzero eigenvalues of pA is bounded by the number 
of zero modes of A^^ particles in N^ fluxes subject to 
the same pseudopotential interaction. For FQH model 
states such as the Laughlin state, the spectrum saturates 
the bound: the number of non zero eigenvalues in pA 
matches exactly the number of quasihole states for Na 
particles and the same number of flux quanta as the orig- 
inal state. The quasihole counting is characteristic of 
each model, thus the PES act s as a fingerp rint. As was 
shown in many previous papera^^EMSIloHlI]^ ^jj^g isotropic 
FQH/FCI state has a low-entanglement energy spectrum 
whose levels saturate the quasihole bound. 

We now ask if the entanglement spectrum of the CDW 
ground-states in the thin torus limit is different than the 
one in the isotropic limit. For the case of our flat-band 
model (and similarly for the FQH thin torus limit pre- 
sented below), we can compute the entanglement spec- 
trum exactly. To make contact with the FQH, we will 
first go to momentum space and then compute the spec- 
trum in that basis (although we will not solve per momen- 
tum quantum number). Let us show as an example, the 
PES for the = 1/3 state in the case when the one-body 
Hamiltonian has n = 1: h^J:i{k) = cos{k)az + sm{k)ax- 
There are 3 CDW states |1), |2) and |3) which are zero 
modes at filling 1/3 of the projected interacting Hamilto- 
nian with on-site and nearest neighbor interaction (which 
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in the projected basis of operators jj becomes a nearest 
neighbor and next nearest neighbor interaction), written 
down in Eq 19 . We can form momentum combinations 



of these states: 



IK) 



(20) 



with K=^K 



f ((A^e + l)mod2) and K being 0,1 or 2. 
The FQH entanglement spectrum is formed by taking the 
superposition of the density matrix of the states at the 3 
K momenta (so that translational invariance is kept). A 
summation over the values of K leads to the identity: 



(21) 



K 



and the density matrix has the same expression in both 
momentum and real space. Tracing out — Na par- 
ticles out of the Nf. electron states, we remain with a 
density matrix for Na particles on the sites of the lat- 
tice. Since the 3 states \j) are orthogonal CDW states 
with electrons fixed at lattice positions, the entanglement 
spectrum breaks up into 3 sectors, and has degenerate 
eigenvalues (this statement is also true for the thin torus 
limit abelian state whose degeneracy is just a center of 
mass degeneracy- for example for the 1/m Laughlin state, 
the spectrum breaks up into m sectors and is completely 
degenerate). The number of states per each sector is just 
the number of ways of choosing Na electrons out of the 
N, 



initial N^ electrons 
values is: 



Na 



The total number of eigen- 



CDW 



Na 



(22) 



All the entanglement energies are exactly degenerate at 
1 /-^cznv^^ • This number is significantly smaller than 
the number of eigenvalues that one would see for the 
Laughlin state, which is 



,Arj ^ 3iVe(3iVe-2A^^-l)! 

NA\{iN,-mA)\ 



(23) 



In the thermodynamic limit, the difference between the 
two counting diverges. The CDW counting is just the 
counting of (1, 3) configurations of Na particles that are 
already present in the CDW groundstate of the system 
while the FQH counting is that of all (1, 3) configurations 
of Na particles in 3iVe orbitals. Hence the entanglement 
spectrum of the CDW does not saturate the bound by a 
large value, and has far fewer levels than the entangle- 
ment spectrum of the FQH state. 

For models without a perfect flat-band limit, the main 
features of the entanglement spectrum remain. The 
CDW-like ground-states will have different entanglement 
spectrum than the topological state. Going from the thin 



torus to the isotropic limit, the remaining levels in the 
entanglement spectrum should come down and mix with 
the CDW-like levels to give the full number of FQH en- 
tanglement levels in Eq 23 . To show this, we plot the 



evolution of the entanglement spectrum of the Kagome 
model from the thin torus limit to the isotropic limit in 
Fig. |3] In the thin torus limit (Fig. [3^), the number 

of entanglement energies is exactly given by J^cdw'"^ 
and the energies are almost all degenerate. The repul- 
sive interaction that we have used is the nearest neigh- 
bor interaction^^. Fig. [sj: and[3^ show how the remaining 
levels mix with the thin torus CDW levels as the aspect 
ratio is made isotropic. Notice that there is a clear en- 
tanglement gap between the low energy physics driven 
by the Laughlin state and higher energy states. Hence 
the entanglement spectrum can be used to differentiate 
between the CDW and the FQH state. A similar calcula- 
tion can be performed on a FQH system at filling factor 
I' = 1/3 on a torus. In that case, we use the threefold 
groundstate of the pseudopotential hamiltoniarP^ that 
generates the Laughlin state. We can then change the 
aspect ratio a — Ly/L^ of the torus to reach the thin 
torus limit and see the effect on the PES. Such results 
are shown in Figs. [3}d, [3ji and [3f. We clearly that both 
the FCI and the FQH PES exhibit the same structure. 
We stress that the PES of the FQH in the thin torus limit 
(almost reached in Fig. [3}d) corresponds to the one of a 
CDW. 



V. THIN TORUS LIMIT IN THE FQH EFFECT 

The physics of the exact fiat-band 1-D models pre- 
sented above is identical to that of the thin-torus limit of 
the FQH. Thus it is highly relevant to look back at this 
case and see which emerging similar structure could be 
look after in FCI. 



A. Hierarchy of bands 

The thin-torus limit of Haldane pseudo- potential 
Hamiltonians is that of electrostatic interactions (or the 
generalized fc-body interaction for the RR state). These 
are just density-density (in the case of Zk RR state, we 
have a fc-body density*^) interactions of density operators 
of the torus momentum j — . . . iV$ — 1. For example, 
for a torus of length , Ly (for simplicity we chose rect- 
angular torus), LxLy — 2'itNip with LLL orbitals in the 
Landau gauge Ay = Bx are indexed by a momentum 
quantum number j: 



1 



' kez 



(24) 



The two-body Haldane 1/3 pseudopotential Hamiltonian 
V{r) = V^(5(r) projected to the LLL, which has the 
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(a) FCI a=l/18 



(c)FCI a=2/9 



(e) FCI a=l /3 



(b) FQHa=l/18 
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(f) FQH a=l/3 
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FIG. 3: PES built from the three low energy states for = 6 
fermions, keeping Na = 3 particles. PES are shown for the 
FCI with iV^ = 18, Ny = 1 (fig. a), ^ 9,Ny ^ 2 (fig. 
c), Nx ~ 6,Ny = 3 (fig. e) and the FQH on a torus with 
aspect ratios a = 1/18 (fig. b), a = 2/9 (fig. d), a = 
1/3 (fig. f). The state counting below the solid green line 
matches the (1,3) Laughlin quasihole counting in Eq[23], the 
state counting below the dashed blue line matches the CDW 
counting in Eq 22 . Notice that the PES of the two orbital 



model in the isotropic case and of the FQH in the thin torus 
limit are identical. 



Laughlin = 1/3 state as the zero energy ground- 
state, can be written using the LLL orbitals in second 
quantized form as Eji,2,3,4 ^ii J2 ja jiCji^CjaCii- In the 



isotropic limit l^^j^j-gj^ can be expressed as an infinite 
sunPS of exponentials, but has large matrix elements 
for large separations ji — j2, etc. It also does not as- 
sume any particularly nice form and cannot be written 
as a density-density interaction, i.e. there is no sense 
in which Vj^j^^j^^j^ with j2 = j3, ji = j4 is larger than 
other matrix elements. In the thin torus limit Ly 
(better described by the ratio a = Ly/L^ — 0, with 
LxLy = constant), a re-organization of the matrix 

elements takes place and several things occur. First, the 
interaction becomes only density-density, and the ma- 
trix elements with j2 ^ J3,ji J4 are exponentially 
suppressed. Amongst the density- density interactions, 
the short-range ones have the highest weight and longer 
range density density interactions decay exponentially. A 
short calculation gives the matrix elements for the Hal- 
dane 1/3 pseudopotential Hamiltonian in terms of the 
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FIG. 4: Left panel : Energy spectrum for FQH on torus and 
the Laughlin f = 1/5 model interaction with — 5 elec- 
trons and aspect ratio a — 0.07. Right panel : Similar energy 
spectrum for the Laughlin v — 1/7 model interaction with 
A^e = 5 electrons and aspect ratio a = 0.095. In both cases 
we observe a zero manifold and several bands associated to 
the (1,7") exclusion principle. The absence of a gap related 
to the (1, 2) for the v = 1/7 case is an artefact of the limited 
immerical accuracy that does not allow to consider smaller as- 
pect ratio without collapsing the (1,7) gap. The plots do not 
show the 5 and 7 fold center of mass momentum degeneracy, 
as the FQH spectrum is exactly degenerate. 



aspect ratio a and number of fluxes iV^ to be: 



V = 



E 



2tt 

aNj, 



8n 



"''■71^71^+2 (25) 



where rii — 



is the density of the j'th orbital mo- 



mentum. This is an exactly solvable model. The exact 
zero modes of this Hamiltonian are the Slater determi- 
nants with occupation number configurations given by 
the rule that no electron can be closer than 3 orbitals to 
another electron (if it is, one or both of the above poten- 
tials are nonzero). This rule is identical to a generalized 
Pauli principle which postulates that there should be no 
more than 1 particle in 3 consecutive orbitals (we call this 
(1,3) Pauli principle). At filling N^/N^ = 1/3, there are 
3 wavefunctions (which are non-interacting Slater deter- 
minants) given by the configurations 100100100100..., 
010010010010 ... and 001001001001. The only difference 
between the zero modes of our flat band model and those 
of the FQH thin torus limit is that they are CDW in real 
space in the flat band model and in momentum space 
in the FQH. What are the excitations of the FQH at 
fillings below 1/3? Notice the separation of scales for 
r ^ in the thin torus Hamiltonian [241 in this limit, 



njnj^2 interaction has a coefficient 4e 



, the 

exponentially 

smaller than that of the UjUj+i interaction. We hence 
expect that the full Hamiltonian has 3 bands of states. 
First the zero modes. Then the first band of excitations 
above the zero-modes of [25] who will be zero modes of 
the HjUj-f-i interaction (which has a large penalty). The 
wavefunctions of these excitations are Slater determi- 
nants made out of configurations which satisfy the Pauli 
principle that there should be no more than 1 particle 
in 2 consecutive orbitals (we call this (1,2) principle). 
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The energy of the excitation spectrum of the band of m = 5, and in Fig. |4b for m = 7. 
states described by this principle starts at 



and ends at {Ng — l);^^e where is the num- 
ber of electrons (we work at filling < 1/3, so the zero 
modes are the FQH quasiholes). For the FQH, the upper 
energy level of the (1,2) excitations is (for any number 
of electrons) smaller than the lower level of the highest 



« 



27T 

aN, 



■ . Hence 



energy band: (iVe-l)^e 

the (1,2) band of excitations (see Fig. |6]d) is separated 
by a gap ^ 

2ir „~ 



from the zero modes and by a gap 



—ij-^ — (Ne — l)-^e °™* from an upper band of 
excitations. The upper band of excitations now is made 
out of all remaining fermionic (Pauli principle (1, 1)) ex- 
citations less the ones in the (1, 2) and (1, 3) bands. The 
Hamiltonian Eq 25 has other terms, including ones that 



cannot be written out as density-density interactions, but 
they are all at least a factor e ""^ smaller than the 
largest term in V. 
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FIG. 5: Energy spectra for A'' = 6 fermions with 18 unit cells 
for different models in the thin torus limit {N^ ~ 1). The 
four models are : the two orbital model with Vi = (a), the 
Haldane model (b), the checkerboard lattice model (c) and 
the ruby lattice model (d). Only the momentum sectors with 
ky < 9 are shown, the other momentum sectors are related 
to these ones by the inversion symmetry. Notice that the two 
orbital models has a highly degenerate manifold for Vi — 0. 



This hierarchy of bands in the thin torus limit is generic 
to the pseudopotentials in FQH, and despite a literature 
search in this extensively researched subject, we were not 
able to find previous mention of it. This statement seems 
to be valid any Laughlin states v = 1/m where we ob- 
serve (1,to) zero modes and then (l,m — 1) up to (1, 1) 
bands, as depicted in Fig. [6]d for m = 3, in Fig. [4^ for 



B. Numerical analysis of FCI models 

We now perform a numerical analysis of the thin torus 
limit of the Kagome model, which, in its isotropic limit, 
has the strongest FCI state. We find a clear separation of 
the states described by the Pauh principles (1,3), (1, 2), 
and (1, 1) in both the case Ne/Ny = 1/3 (see Fig(6]a). A 
similar result is also observe for N^/Ny < 1/3, i.e. when 
quasiholes are considered (see Fig[7]a). 

We now conjecture that good FCI states have, in the 
thin torus limit, a similar separation of scales as the FQH. 
The FCI projected Hamiltonian is similar in spirit to the 
FQH hamiltonian: when projected, in the isotropic limit, 
it cannot be written as density-density interactions of any 
range. In the thin torus limit, at least for models where 
the band is flat without the need for long-range pro- 
jectors, the FCI projected Hamiltonian becomes short- 
ranged and writable in terms of just density-density in- 
teractions For FCI Hamiltonians which when projected 
in the thin torus limit, exhibit a separation of scales sim- 
ilar to that in the FQH we observe strong FCI states in 
the isotropic limit. This is the criterion which we conjec- 
ture is essential in finding FCI states on the lattice, and 
we support it below with numerical data. Fig. [5] shows 
the energy s pectrum for several models. The Rub> *^^ l ^^ l, 
checkerboar d^ l ^ l ^ ^^^, and Haldane modeP^, which ex- 
hibit decreasingly strong FQH states in the isotropic limit 
for fermions, have a decreasingly clear separation of scales 
in the thin torus limit. The ruby model, with one of the 
strongest FQH Laughlin state^^ has a clear manifold of 3 
zero modes (only 2 are shown, the third one being the in- 
version symmetric of the one ai ky = 3, and shows up at 
ky = 15), then a manifold of (1, 2) Pauli principle states 
clearly separated from both the zero modes and from the 
highest energy manifold of (1, 1) states. The bandwidth 
of each of these manifolds is small, and they are clearly 
defined. The Haldane model still shows a similar separa- 
tion of scales, but the bandwidth of each manifold is now 
large and unlike the FQH. The two-orbital model with 
just Vq does not show such a separation (for example, 
there is no difference between the (1, 3) and (1, 2) modes 
which appear at zero energy), as we can prove from an- 
alytics (we need the Vi to lift the degeneracy of (1, 2) vs 
(1,3) modes), and the isotropic state in the two orbital 
model is not a Laughlin state at this lattice size. 



VI. CONCLUSION 

In this paper we have analyzed the thin torus limit 
of several models which are supposed to exhibit a FCI 
state. The thin torus limit of lattice models was de- 
fined as the Nx — 1 limit of an x Ny site lattice 
of electrons in the presence of any Hubbard-like interac- 
tion. We found a series of one-body models with finite- 
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FIG. 6: Energy spectra for A'^ = 6 fermions for the Kagome 
lattice with 18 unit cells (upper panel) and FQH on the torus 
geometry with N<s> = 18 flux quanta and aspect ratio 1/18 
(lower panel). In the FQH case, we only display the mo- 
mentum sectors that cannot be related by the center of mass 
degeneracy. The small dispersion of the low energy manifold 
on the FQH spectrum is due to the numerical accuracy. 



range hoppings which for isotropic aspect ratio exhibit 
nonzero Chern numbers in the nontrivial side and which, 
in the thin-torus limit, exhibit perfectly flat bands. In 
the thin torus limit, the nontrivial side of the insulator 
exhibits dimerization, trimerization, or more generally 
n-merization etc. Our thin torus limit one-body Hamil- 
tonians can also be thought as inversion symmetric one- 
dimensional insulators with half charge polarization on 
the edge. In our toy models, due to the short-range form 
of the projection operator, we can diagonalize exactly 
Hubbard-type interactions in the lowest band. We show 
that projecting the interaction to the lowest band in- 
volves an increase in its range by exactly the distance 
of the n-merized bond in our one-body model. In our 
toy models the interactions can be tuned so that when 
projected they give exactly the correct ground-state de- 
generacy for the state we are considering. We choose the 
interactions so that the spectrum matches qualitatively 
that of the thin-torus limit of the FQH. Due to their sim- 
plicity, we can exactly solve the full spectrum of such in- 
teractions, and show that it separates into bands of states 
spanned by Slater determinants in real space satisfying 
a certain exclusion statistics hierarchy. For example, for 
the interaction responsible for the 1/3 filling Laughlin 



FIG. 7: Energy spectra for TV = 6 fermions for the Kagome 
lattice with 20 unit cells (upper panel) and FQH on the torus 
geometry with N<s> = 20 flux quanta and aspect ratio 1/20 
(lower panel). This situation corresponds to the Laughlin 
state with two added quasihole. In the FQH case, we only 
display the momentum sectors that cannot be related by the 
center of mass degeneracy. The small dispersion of the low 
energy manifold on the FQH spectrum is due to the numerical 
accuracy. 



state, in the thin torus limit the spectrum is made out 
of 3 bands: the zero modes are exactly the ground-state 
and quasiholes of the Laughlin state satisfying a Pauli 
(1,3) principle that disallows more than 1 particle in 3 
consecutive orbitals. The middle band satisfies a (1,2) 
Pauli principle, while the highest energy band satisfies a 
(1,1) (usual fermionic) Pauli principle. In the thin torus 
limit, all these states are CDW (or its excitations), which 
shows pure counting of ground-states and excitations is 
not enough to fully determine the topological nature of 
a state. We then compute analytically the entanglement 
spectrum of the degenerate manifold of states and show 
that the entanglement spectrum of the CDW state and 
that of the isotropic FCI state differ considerably, both 
qualitatively and quantitatively, and can be used to dif- 
ferentiate between the states. We identify the CDW and 
FCI structures in the entanglement spectrum and show 
how the CDW spectrum evolves into that of the FCI upon 
going from the thin-torus to the isotropic limit. We then 
numerically analyze several FCI models and present nu- 
merical evidence for the conjecture that strong isotropic 
FCI states appear when there exists a separation of scales 
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in the thin torus hmit similar to that of the FQH. 
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